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Fast three-dimensional (3D) imaging requires parallel optical slicing of a specimen with an efficient 
detection scheme. The generation of multiple localized dot-like excitation structures solves the problem of 
simultaneous slicing multiple specimen layers, but an efficient detection scheme is necessary. Confocal theta 
detection (detection at 90° to the optical axis) provides a suitable detection platform that is capable of 
cross- talk-free fluorescence detection from each nanodot (axial dimension —150 nm). Additionally, this 
technique has the unique feature of imaging a specimen at a large working distance with super- resolution 
capabilities. Polarization studies show distinct field structures for fixed and fluid samples, indicating a 
non-negligible field-dipole interaction. The realization of the proposed imaging technique will advance and 
diversify multiphoton fluorescence microscopy for numerous applications in nanobioimaging and optical 
engineering. 

Optical sectioning 1 and point-by-point scanning are required for most existing fluorescence microscopy 
techniques. This technique have the advantage of being non- destructive, and they offer the flexibility of 
slicing along the desired sample layer. Confocal laser scanning microscopy and multiphoton fluor- 
escence microscopy are routinely used to make optical scans of a sample and stack them for 3D reconstruction 
of biological structures 2 ' 3 ' 4 ' 5 . Selective plane illumination microscopy (SPIM) is a similar technique in which single 
planes are illuminated sequentially to obtain sectional images 7 . Recently, there has been growing interest in 
multiphoton super-resolution scanning techniques 8,910 . Despite the rapid advances in light microscopy, sim- 
ultaneous visualization of multiple layers (that includes both excitation and detection) has not yet been achieved. 
Such a capability would have numerous applications, including, fast imaging (by parallelizing data acquisition), 
and data correlation across multiple layers of the specimen. For example, simultaneous data from the nucleus and 
the cell membrane could be recorded and correlated to extract vital information about the internal function of a 
cell. Furthermore, existing fluorescence imaging techniques offer limited imaging depth-, and poor axial reso- 
lution, and they suffer from photobleaching 11 ' 12 ' 13 ' 14 ' 15 . Proposed technique holds the promise of reducing the 
photobleaching that is prevalent in biological imaging and offers better depth-imaging capability. For this reason, 
we have previously proposed an excitation technique for obtaining multiple localized field patterns using an 
optical mask 1617 . The proposed excitation scheme is straight forward, but the detection scheme is tedious, 
especially when distinguishing the emission from the nanoscopic dots while simultaneously avoiding cross-talks 
between them. The goal is to excite multiple specimen layers that are few hundred nanometers apart, and to 
selectively detect fluorescence from the nanodots almost simultaneously. 

In this report, we propose an orthogonal detection scheme for multiphoton multiple excitation spot optical (M- 
MESO) microscopy, in which the emission is collected by an objective that is placed perpendicular to the optical 
axis. This orthogonal detection was initially developed by Voie et al. 1 and further developed by Stelzer et al. 18 . It 
should be noted that, a long working distance objective lens has a low numerical aperture and the converse is also 
true. For the proposed M-MESO microscopy, confocal theta detection offers an effective solution because this 
configuration combines the advantage of high resolution and long working distance. Here, we show that, the 
emission from a single nanodot can be effectively collected without introducing cross -talks from its neighboring 
nanodots. Orthogonal detection has the following additional features compared to conventional confocal detec- 
tion: (1) the system axial resolution is effectively the lateral resolution of the detection system, (2) most of the 
excitation light is cut-off because of the perpendicular detection (and the background light is substantially 
reduced as a result), and (3) the detection arm is independent of the excitation arm, which provides additional 
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flexibility for fluorescence detection without altering the excitation 
mechanism. We show unique (cross-talk free) identification of an 
individual nanodot from a series of nanodots by utilizing a confocal 
theta detection scheme. In principle, simultaneous detection from all 
of the nanodots can be obtained by introducing a distorted diffrac- 
tion grating in the detection path of the imaging system 19 , followed by 
the charge- coupled device (CCD) detector. Interestingly, this detec- 
tion technique reports nano-particle tracking of approximately 
lOnm along the z-axis 19 . Alternately, a 50/50 beam-splitter with an 
appropriately positioned pinhole can be used in the detection path- 
way 20 . The proposed imaging technique is a boon to multiphoton 
microscopy, in which the poor axial resolution of the multiphoton 
imaging process limits the advantage of thin optical slicing. 

Results 

The conventional laser- scanning imaging technique is limited by the 
classical resolution limit and suffers from slow data acquisition as a 
result of the point-by-point scan. These shortcomings encourage the 
development of fast super- resolution fluorescence imaging tech- 
niques. Ideally, the technique should have simultaneous multi-layer 
excitation and cross-talk-free detection for diffraction-unlimited 
multilayer imaging. These attributes enable several applications, 
including fast imaging and multi-layer crosscorrelation studies. 
Here, we study a new imaging system (see Figure 1) and specifically 
characterize the system point spread function (PSF) under the iso- 
tropic emission model with randomly polarized light illumination. 
Illumination, detection and the system PSF for the proposed M - 
MESO system are shown in Figure 2. The axial resolution of the 
proposed system is improved by a factor of approximately 4 over 
the classical resolution limit, while there is no improvement in the 
lateral resolution because the perpendicular detection occurs along 
the lateral axis of the illumination objective. As a result, the axial 
resolution of the total imaging system is essentially determined by the 
lateral resolution of the objective 0 3 . Thus, the multi-dot illumina- 
tion PSF (km) multiplied by the detection PSF (h det ) gives the system 
PSF, which is the desired single-dot (see Figure 2). The correspond- 
ing axial and lateral resolutions are approximately, 210 nm and 
150 nm, respectively (see Figure 2). Moreover, a system PSF volume 
of approximately 27.7 X 10~ 3 femto - lit is obtained. By scanning the 
objective 0 3 , the fluorescence from all of the nanodots can be 
recorded, in principle. Moreover, we have carried out studies on 
the single photon variant of the proposed MESO imaging system. 
This study facilitates the direct comparison of the single- and multi- 
photon MESO imaging techniques. The excitation, detection and 
overall system PSF are shown in Figure 3 for otm = 45° and ot det = 
45°. Notably, the density of the nanodots doubles because of the 
smaller wavelength employed for the single-photon MESO system. 
The size of the individual nanodot (determined by its full -width 
half-maximum) is small (lateral resolution r xy = 180 nm and axial 
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Figure 2 | M-MESO system PSF for a m =45° and a det =30°. 

resolution Az = 120 nm) for single photon MESO compared to 
the M-MESO system, for which r xy = 210 nm and Az = 150 nm. 
Increased crosstalks have been observed in the single photon variant 
compared to the multiphoton variant. This increase is evident in 
Figure 4, which shows the intensity plots obtained along predefined 
lines passing through the system PSF. These plots indicate prominent 
sidelobes for single-photon MESO compared to M-MESO, which 
can be attributed to the doubling of the nanodot density as a result 
of the shorter excitation wavelength used for single photon MESO. 
The shorter wavelength modifies the interference condition, thereby 
altering the distance between nearby nanodots. As a result, the detec- 
tion PSF picks up fluorescence from the nearby nanodots (both 
sides). This difference clearly illustrates the advantage of M-MESO 
over its single photon variant, in addition to the other advantages of 
the multiphoton excitation process. 

In a high resolution multiphoton imaging system, the advantage of 
axial super- resolution at long working distances is of paramount 
importance. For this reason, we have investigated the system 
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Figure 1 | Optical configuration of the proposed M-MESO microscopy. Figure 3 | Single photon MESO PSF for a m = 45° and a det = 45°. 
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Figure 4 | Intensity plots along the X Z-plane. The inset shows the 
corresponding single- (1PE) and two-(2PE) photon MESO PSFs. 

performance at a small illumination aperture angle (am = 18°) and 
varying detection aperture angles (oc det ). At low numerical aperature 
of the detection objective (oc det = 30°), a sidelobe appears because of 
the increased crosstalks from both of the neighboring excitation dots 
(indicated by red arrows in Figure 5(a)). Further evidence of this 
cross-talk is shown in the intensity plots (see Figure 5(b)). The effect 
is substantially minimized at oc det = 60° (see Figure 5(b)), and com- 
plete elimination of sidelobes occurs at high NA detection (ot det = 
72°). This trend occurs because of the high lateral resolution of the 
detection system. For the proposed imaging system, the axial reso- 
lution is determined by the product of the individual excitation 
nanodot and the lateral PSF of the detection objective, 0 3 . Overall, 
detection at large aperture angles provides information from indi- 
vidual excitation dots. Because the high NA illumination objective is 
practically helpful, we have also studied a different configuration of 
M-MESO system, in which we varied the illumination aperture angle 
(oc m = 30°, 60°, 72°) and fixed the detection aperture angle oc det at 45°, 
as shown in Figure 6. The nanodots are observed to elongate along 
the optical axis for large dm, thereby increasing the distance between 
two nearby nanodots. As, a result, the crosstalk from the neighboring 
nanodots is completely eliminated (see the intensity plots in 
Figure 6). This configuration is beneficial for artifact-free high reso- 
lution imaging. To elucidate the multilayer imaging capability of the 



M-MESO system, we have detected nanodots at varying depths (z- 
layers). Figure 7, shows the detection of target nanodot situated at 
depths of 540 nm and 1.56 /urn, along with the detection of a central 
nanodot (at z = 0 nm). The target nanodots are indicated by red 
arrows in the excitation PSF (h exc ). Simultaneous detection from all 
nanodots can be obtained by using appropriate optical elements, 
such as, a diffraction grating in the detection path of the imaging 
system 19 . Such a capability is useful for near-simultaneous indepth 
multilayer imaging of biological specimens. 

The optical transfer function (OTF) of an imaging system provides 
bandwidth information. The OTF is obtained by the convolution of 
the excitation and detection OTFs. The OTFs for both the single- 
and multi- photon versions of the proposed imaging system are 
shown in Figure 8. The OTF peaks at zero frequency (k xy = 0, k z 
= 0). Missing bands were observed along the optical axis for single 
photon MESO, while the M-MESO imaging system shows continu- 
ously filled bands. The role of multiphoton excitation is to fill the 
missing bands in the three dimensional OTF. Gaps in the OTF 
indicate missing object spectrum information. Overall, better-filled 
bandwidth improves the efficiency of 3D image reconstruction 
because gaps in the OTF pose problems to the ill-posed image recon- 
struction process. Additionally, the two distinct sidelobes in the 
single photon MESO indicate substantial crosstalk that is missing 
in M-MESO, thereby illustrating the advantage of the M-MESO 
technique over its single photon variant. 

Practical situations demand knowledge of the field distribution for 
both fixed and fluid samples. In fluid samples, the dipoles are free to 
rotate whereas the dipoles are static in fixed samples. As a result, the 
same field structure cannot be assumed for both cases, and the dipole 
orientation must be taken into account for field calculations. The 
interaction between the electric field and the molecular dipoles must 
also be considered. Scalar-diffraction theory has proven to be valu- 
able for understanding the structure of the electric field of various 
optical components, such as lenses. However, the theory is not 
applicable for high numerical aperture (NA) objective lenses, which 
requires the tight focussing of light for fluorescence-based applica- 
tions. Moreover, the theory lacks vectorial features of the electric 
field, such as, polarization, and it is unable to explain polarization 
dependent effects. This scenario requires the complete clasical theory 
of light, i.e, the vectorial -diffraction theory. Vectorial theory has been 
employed to quantify the electromagnetic field structure arising from 
both field polarization and dipole orientation for the proposed 
M-MESO imaging system. Figure 9 (top) shows the lateral (XT) field 
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Figure 5 | (a) System PSF for varying detection aperture angle a det = (30°, 60°, 72°) and illumination aperture angle a m = 18°, (b) Corresponding 
intensity plots. 
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Figure 6 | (a) System PSF for varying illumination aperture angle a m = (30°, 60°, 72°) and fixed detection aperture angle a det =45° 
(b) Corresponding intensity plots. 
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Figure 7 | Nanodot detection at varying z-layers: 0 nm, 540 nm and 1560 nm. Red arrow indicates a target nanodot in the excitation PSF. Both the 
x- and z- axes extends from — 1.92 fim to 1.92 fim. 
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Figure 8 | Optical transfer functions for both the single photon MESO and M-MESO imaging systems. Red arrows indicates missing bands in 3D OTF 
of single photon MESO system. 



distribution for free-rotating dipoles illuminated by both random- 
polarized (RP) and linear-polarized (LP) light. The region of uniform 
excitation is large for LP light, but for RP light illumination, the 
uniformity is split into three prominent parts (see Figure 9, top-left). 
This result indicates that, LP light is beneficial for applications 
requiring uniform excitation. Figure 9 (top row) and (bottom row) 
show the field structure of a free rotating (representing fluidic sam- 
ple) and x- oriented dipole (representing fixed sample), respectively. 
Irrespective of the polarization of the electric field, the field structure 
is significantly different and depend strongly on the dipole orienta- 
tion, thereby eliciting the prominance of the field- dipole interaction. 
Understanding the effect of light polarization on the orientation of 
dipoles is essential for the realization of a precision super-resolution 
imaging system. 

Finally, the proposed system (M - MESO) is compared with two- 
photon excitation (2PE) and two -photon excited 4pi (2PE - 4pi) 



imaging systems. Confocal detection perpendicular to the optical 
axis is employed for all cases. The illumination (hm), detection (hd et ), 
and the system PSF (h sys ) are shown in Figure 10 for all imaging 
modalities. High resolution is evident in the 2PE system because of 
its confocal theta detection. Comparatively, the resolution is much 
better for 2PE - 4pi and M - MESO. However, the 2PE - 4pi system 
can generate only one super-resolution dot, whereas, the proposed 
M - MESO system can generate several such high-resolution dots 
of dimension Az ~ 150 nm. Additionally, the M-MESO system is 
capable of fluorescence detection from each individual nanodot 
without any cross-talks from neighboring nanodots. 

Discussion 

In principle, it is possible to further reduce the layer thickness with 
the use of a high NA objective. However, this alteration reduces the 
working distance thereby making it difficult to realize in practice. On 
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Figure 9 | Field structure of four different combinations of field polarization and dipole orientation: (top) RP and LP field on freely rotating dipole, 
(bottom) RP and LP field on x-oriented fixed dipoles. 



the other hand, it is also possible to increase the number of excitation 
dots by decreasing the NA of the excitation objective, but this vari- 
ation enhances the cross-talks between the nearby nanodots, 
destroying the unique detection of fluorescence from each nanodot. 
Consequently, a balance must maintained for practically realizing the 
proposed imaging technique. Nonetheless, the parallel excitation of 
multiple specimen layers using nanodots is exciting. A spatial filter 
based multidot illumination, coupled with an orthogonal detection 
PSF, shows considerable tolerance for crosstalk-free fluorescence 
collection. This excitation technique is promising because the spatial 
filter pattern determines the shape and size of the system PSF. In 
addition, such filters enable diffraction-unlimited lateral resolution, 
sub -femto -liter volume excitation and the generation of Bessel-like 
beams 28,29 . As a result of the convenience and new possibilities assoc- 
iated with M-MESO, we anticipate that this technique may become a 
powerful tool for numerous applications in Biology, nanoscale 
imaging and optical engineering. These applications may include 
multiple-trap optical tweezers and multi-dot fluorescence correla- 
tion spectroscopy. 

Scattering and anisotropy have a significant effect on the quality of 
the image and on multi-dimensional image reconstruction. These 
effects are minimized for two-photon imaging that uses, near infra- 
red radiation (IR) for multiphoton excitation because the use of IR 
minimizes scattering 22,23,24 . The quadratic dependence on the excita- 
tion light intensity for two -photon excitation makes scattering less 
important because the excited fluorophore is spatially localized in the 
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Figure 10 | Comparison of 2PE, 2PE-4pi and M-MESO PSF at an 
illumination and detection aperture angle of 45 ° . 



highest intensity regions of the beam. This mechanism, which is 
associated with two -photon excitation, leads to lower background 
and less photobleaching. In Rayleigh scattering, which is the most 
prominent scattering process for particles much smaller than the 
wavelength of light, the amount of scattered light scales inversely 
as the fourth power of the wavelength. Therefore, at an excitation 
wavelength of 910 nm (used for 2PE), scattering intensity is approxi- 
mately 8 times less than for 488 nm (used for 1PE) excitation. This 
reduction in scattering increases the depth of penetration in thick 
biological specimens. Biological specimens are non-uniform and 
have variable indices of refraction, which causes multiple scattering 
of light in various directions. Specifically, in fluorescence imaging, 
the incident light is scattered to varying degrees before it reaches the 
focal plane. The resulting fluorescent light also undergoes similar 
scattering as it travels through the specimen toward the detector. 
Both of these scattering effects reduce the collected fluorescence 
signal, resulting in numerous image artifacts and hampering the 
image resolution. On the other hand, for particles a few microns in 
size, Mie scattering is dominant, and the scattering efficiency 
(defined in terms of the scattering cross-section and anisotropy) of 
the specimen determines the image resolution 25 . A decrease in the 
scattering efficiency implies a reduction in the number of multiple 
scattering events, resulting in increased signal collection by the 
detector and high image resolution. TPE is reported to have an image 
resolution that is two orders of magnitude higher than that under 
1PE 25 . 

Methods 

M-MESO excitation is employed for the simultaneous multi-layer slicing of biological 
specimens. An efficient detection scheme is used to collect photons from the target 
nanodots to realize the complete imaging system (see Figure 1). We recognize con- 
focal theta detection as the most suitable detection scheme for the proposed imaging 
system. Phase-matched light is structured by the spatial filter (imprinted on the 
objective lenses O x and 0 2 ) and focused to the common geometrical focus of a 4n- 
geometry 9 ' 26 ' 27 (47r-geometry is extensively used in fluorescence microscopy for axial 
super-resolution 9 ). This focusing results in the formation of a series of bright nano- 
dots because of interference along the optical axis. Aperture engineering on a single 
lens configuration with the proposed spatial filter produces a Bessel-like depth- of- 
focus PSF 28 ' 29 . Other spatial filters are also proposed for high resolution imaging 30,31 ' 32 . 
Fluorescence from these tiny nanodots is collected by the objective 0 3 , which is 
positioned orthogonally to the optical axis (see, Fig 1). For the proposed study, the 
two-photon excitation and emission wavelengths are chosen to be 910 nm and 
520 nm respectively. The excitation and emission/detection wavelength for the single 
photon variant are chosen to be 488 nm and 520 nm, respectively. The 3D image 
consists of 128 X 128 pixels and 128 slices, with the geometrical focus located at z = 0. 
A sampling rate (well above the Nyquist sampling) of about 30nm is chosen along 
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both the axial and lateral axes. To generate excitation nanodots, a transmission 
window of Ad = (cc-6{) = 5° is chosen for all cases (see Figure 1). 

A schematic diagram of the proposed imaging system is shown in Figure 1 . 
In the Cartesian coordinate system, the x-,y- and z- components of the electric field 
for randomly polarized light illumination are given by 21 , [E x , E yy E z ] = [ — iA(I 0 +I 2 
cos(20)), — iAI 2 sin(20), — 2Al x sin(^)], where A is the proportionality constant 
representing the amplitude of the incident electric field. Assuming substantial 
interaction between the incident electric field and the random-oriented 
dipoles, the excitation PSF of the M-MESO system is, hm(x,y, z) = 

\^x,y,z)+^x,y, -z)\ 4 =\\Rel 0 \ 2 + 2\Reh\ 2 + \ReI 2 \ 2 ] , where I Q12 = 

L) 

Seto 5 (*)Go,i,2(*)cos 1 / 2 9 e V sin2 w dO are modified diffraction integrals with, 
[Go(*), GiH, G 2 (*)] = [sin 6(1 + cos 0)/ o (w), sin 2 6h(w), sin 6(1 - cos 6)} 2 (w)], 
where, w = (v sin 6/sina). The function, = H[6— 6 X ] — H[6— a] represents the 
spatial filter with cutoff angle 6 U for which the trans mission window is A6 = 5°; H is 
the Heaviside function. u=^f z sin 2 a and v = ^/x 2 +y 2 sin a are the longitudinal 
and transverse optical coordinates, respectively 21 , and am is the illumination semi- 
aperture angle. The numerical aperture of the objective lens in air is defined as, sin a, 
where a is the semi-aperture angle. For single photon MESO, the illumination PSF is, 
km = \ReIo | 2 + 2\ReIi | 2 + \ReI 2 1 2 . We employ a confocal theta detection scheme, in 
which the detection is carried out in the orthogonal plane, as shown in Figure 1 . The 
orthogonal plane of detection is represented by the following transformation: (x' = 
—z,y' = y, z' = x). This detection scheme has the advantage of high resolution at long 
working distance. The isotropic emission model is assumed with randomly polarized 
light excitation. The corresponding detection electric field is, hd e t(x',y',z') = 
\E X ' | + \Ey I + \E Z > | , where the integration J 0 ,i,2 = Sq=o (' ' on the aperture-free 
objective (0 3 ) is carried over the detection semi-aperture angle a det . Overall, the 
system PSF of the proposed imaging system is given by, h sys (x, y, z) = h m (x, y, z) X 
h det (-z,y,x). 

The incident electric field gives rise to an oscillating dipole (p), which in turn 
produces a far-field electric field given by, E = — r x (rxp). Assuming paraxial optics 
and using the two-dimensional (2D) Fourier transform, one can obtain explicit 
expressions for the calculation of the time-averaged electric energy density, following 
the approach of Sheppard et al 33 . Here, the field is calculated for four different cases. 
Initially, the RP light on a freely rotating dipole is considered, for which the field 
distribution is, h RP j ree rot = |j Rel 0 1 2 + 2 1 Rel\ | 2 + | Rel 2 1 2 j . For x-oriented dipoles, 
the field distribution becomes h^p^- or i ent = ^\ReI 0 1 2 + 2\Reli | \ReI 2 |cos(20) + 
\ReI 2 1 2 ] 2 . Considering LP light on a randomly oriented dipole, the field distribution is 
given by, h LP j ree - rot = [\ReI 0 \ 2 + 4\ReI l | 2 cos 2 ^ + \ReI 0 | 2 cos(20) + \ReI 2 | 2 j , and 
the field distribution becomes, h L p^- or i ent = ^\ReIo | 2 + 2\ReI 0 1 \ReI 2 |cos 2 ^ + 
\ReI 2 1 cos 2(j)] 2 for a fixed dipole orientation along x. 

In the case of fast imaging, to scan M slices of an N X N image, a laser scanning 
system requires MXN 2 scan-points and a total time of (M X N 2 ) At. The generation of 
an excitation PSF with M axial dots, allows in-principle the data acquisition rate to be 
increased by a factor of M, thereby reducing the total acquisition time by a factor of M 
At. The proposed technique can be useful for fast imaging of biological specimens 
with reduced background and photobleaching effects. 
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